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Abstract 

We propose an analog of the familiar gap equation for the case of four-quark conden- 
sates at finite temperature. The condensates of interest correspond to scalar, vector, 
psudoscalar, axial vector, and tensor Dirac structures. Working with correlators at 
zero chemical potential without factorization, we arrive at coupled equations for these 
four-quark condensates and the masses of certain light-quark mesons. We study the 
temperature dependence of the four-quark condensates and masses; in one of our mod- 
els, factorization of the four-quark condensates is shown to be increasingly violated as 
the temperature is increased toward Tc- The 2"'""'" tensor mesons a2(1320)-f2(1270) are 
identified as especially sensitive probes of the four-quark condensates. 
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1 Introduction 



The early universe is believed to have evolved through a stage beginning about 
100 ps and ending about one /x sec after the big bang, during which time leptons had 
acquired their masses and strongly interacting matter existed primarily as a plasma of 
quarks and gluons. As the universe further expanded and cooled, when the temperature 
reached the value T = Tc ~ 200 MeV matter began to undergo a phase transition 
as quark condensates and possibly hadrons began to form. Upon further expansion 
and cooling, this transition ran to completion with all quarks and gluons then having 
combined into hadrons. 

Our paper concerns the physics that developed as the universe passed from T c^Tc 
to r = 0. The relevant issues are identified by focusing on mesons in such a medium. A 
large body of work on vector mesons at finite T already exists, having been motivated 
by the fact that vector mesons can couple directly to photons providing an experimental 
lepton-pair signal of the masses and widths of these particles. However, without the 
factorization assumptions upon which most of this work relies, and which we describe 
in detail below, all the mesons become coupled through four-quark condensates. In 
the present paper we describe this physics without factorization, deriving a set of 
coupled equations for the meson masses and for the relavant four-quark condensates, 
and giving solutions with various scenarios for the finite-T spectral densities and two- 
and four-quark condensates. 

The quantitative characterization of this sequence of events lies beyond the ca- 
pabilities of theoretical analysis at present. However, because of its high interest, 
experimental programs have been planned to create similar conditions in the labora- 
tory. It is the explicit aim for the RHIC (Relativistic Heavy Ion Collider) project at 
Brookhaven National Laboratory (BNL) to initiate a transition to matter at T pa 200 
MeV or greater, the quark-gluon plasma, within a large enough volume, and for a 
sufficiently extended time, to allow the study of its properties. This is envisioned to 
occur when the matter of large nuclei is compressed and heated as a result of violent. 
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head-on collisions between nuclei that have been accelerated to relativistic energies and 
allowed to collide in the RHIC storage ring. Encouraging evidence from preliminary 
experiments at the AGS at BNL and the SPS at CERN that a hot, dense fireball is 
formed in relativistic heavy ion collisions [l| keeps interest high in the RHIC program. 
The experimental determination of the nature of this transition is expected to guide 
the development of theoretical methods to treat the same physics. 

In this paper we will develop analytical tools capable of addressing the temperature 
dependence of the four-quark condensates characterizing the quark-gluon plasma that 
may be created in this fashion. We are particularly interested in the four-quark conden- 
sates because they are fundamental quantities that characterize the nonperturbative 
vacuum in quantum chromodynamics (QCD) and are needed within the framework of 
QCD sum rules [Q] to describe properties of both mesons (for a review, see Ref. [^]) 
and baryons Q. They are especially important in the case of the light-quark mesons, 
which are the lowest excitations of the physical vacuum. The seminal work on this 
problem at T 7^ was done by Bochkarev and Shaposhnikov, Ref. [Q. They investi- 
gated an important four-quark condensate, namely that contributing to the mass of 
the p meson, using QCD sum rules extended to finite temperature. 

Bochkarev and Shaposhnikov introduced a dynamical model for the temperature 
dependence of this four-quark condensate, thus avoiding the so-called factorization 
approximation, according to which each four-quark condensate becomes trivially pro- 
portional to the square of the familiar quark condensate, assuming saturation by the 
vacuum. A number of other authors have also calculated the mass of the rho meson 
using QCD sum rules at finite T |^, ^, generally using factorization. For the related 
physics of hadrons at T = but finite density it has also been seen § that the fac- 
torization of the four-quark condensates is a major question, and that, e.g., in nuclear 
matter the theoretical result for the nucleon mass is criticaly dependent on the fac- 
torizaton assumption. In a study of the A(1236) in nuclear matter it has been shown 
that factorization of four-quark condensates is not valid for hadrons in a medium of 
non-zero density |p. At zero baryon density there is evidence that it is a bad approxi- 
mation ||l^; however, because it is simple and convenient the approximation is widely 
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used. To our knowledge, the present work is the first attention given to the study of 
four-quark condensate in dynamical models since the original work of Bochkarev and 
Shaoposhnikov. 

An important set of four-quark condensates are those of the form 

= (0| : qWqqreq : |0) (1) 

where are the SU(3) Gell-Mann color matrices and T corresponds to the set of five 
independent Dirac structures T = 1, 75, 'J^,^5j^, and a^^, with an implied contraction 
over the Dirac indices. It is particularly relevant for the present work that these con- 
densates make major contributions to the masses of the light-quark mesons generated 
by the corresponding five quark currents, 

= q{x)Tq{x), (2) 

where q(x) are the u,d-quark fields contracted over their color labels. The mesons 
generated by a particular current J^{x) are those that are excited from the vacuum by 
the action of the current, i.e., those that satisfy 

< 0\J^ {0)\meson(T) > = ^r^'rCp) / 0, (3) 

where ^'r(p. A) is a wave function characteristic of the meson and gr is a structure 
parameter. We will often refer to the mesons and other quantities associated with the 
five different Dirac T by s, ps, v, av, and to-, respectively. 

In our discussion of the mesons, we will also be interested in the axial, tensor 
and axial tensor mesons because their masses depend on the same set of five four- 
quark condensates. The currents required for these mesons |^ are somewhat more 
complicated than the forms given in Eq. (^). One should also note that to determine 
the mass in the scalar (0++) case, the current must include scalar glueball mixing; this 
is discussed in detail below. The most important mesons associated with each current 
will be specified below. 

2 QCD Sum Rules for Mesons 
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The use of QCD sum rules for our purpose is suggested because this framework 
makes a relatively simple and explicit connection between the four-quark condensates 
and light meson masses. QCD sum rules have a number of drawbacks, however, such 
as uncertainty in how to treat the continuum, which becomes more and more serious 
issue as the temperature is raised. For this and other reasons, we do not aim to make 
predictions of QCD but rather develop QCD sum-rule based models, in which param- 
eters are introduced to account for the temperature dependence of the continuum. We 
find our results reasonable, and we believe our ad hoc assumptions are capable of being 
improved as experimental results warrent it. 



2.1 General Considerations 



In the method of QCD Sum Rules II , the propagator of a hadron is represented by 
a two-point function, usually called a correlator, defined as 

lf{p) = i J d^xe"-P <0\T[J^{x)J^{0)]\0>, (4) 

where for light-quark mesons the currents J^(x) are given in Eq. (^) and other forms 
specified below. This differs from the two-point function of standard field theories 
in that the J^(x)'s are composite. Similar correlators are introduced in the study of 
baryons, but in the present work we focus only on the mesons for reasons discussed 
above. 

A QCD sum rule for a quantity of interest is possible when the correlator is saturated 
by hadronic states, i.e. by using Eq. (^) in Eq. (Q), producing a dispersion relation in 
which the quantity of interest appears explicitly. The expression resulting from this 
method of evaluation is often called the right-hand side (R.H.S.) of the sum rule. One 
form for the dispersion relation, used in the present work, is: 

uWsip) = r ^^^^ (5) 

with 

Xa 
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(a) perturbative 
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(b) quark condensate 
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(c) mixed condensate 



Figure 1: Quark propagator. 

where \maip,a)) of 4-momentum {Ea,p) is a complete set of intermediate states of 
quantum numbers A, including the continuum. 

To make a connection to QCD, the correlator is then evaluated microscopically by 
expanding in operators of increasing dimension using an operator product expansion 
(o.p.e). The resulting evaluation is often called the left-hand side (L.H.S.). For the 
quark propagator in space-time, 

Sg{x) = <0\T[q{x)m]\0>, 

= 5f^-(x) + 5f^-(x), (7) 

with 

Sf^-(x) = < 0| : g(x)g(0) : |0 >, (8) 

where S^''^'{x) is the usual perturbative quark propagator, illustrated in Fig. la, and 
S^'^-{x) is the nonperturbative quark propagator, which vanishes in the perturbative 
vacuum. The operator product expansion gives 

= ^ <0| igg: |0>+^ <0| |0>+... (9) 

in terms of the quark condensate, whose value is given by 

(0| : qq : |0) = -(0.25 GeV)^ (10) 

a mixed condensate, and so forth as illustrated in Figs, lb, Ic. The justification for 
the o.p.e. is that the evaluation of Il^{p) is done using a Borel transformation with 
the Borel mass in the range of about 1 GeV, allowing a short-distance expansion. 

The microscopic derivation of the correlator Il^{p) is completed by using the ex- 
pansion of Sg(x) given in Eqs. (0,^), leading to the perturbative term illustrated in 
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(a) perturbative (b) gluon condensate (c) four-quark condensate 



Figure 2: O.P.E. for meson correlator. 



Fig. 2a. To this is added terms with gluons attached to two or more quark propaga- 
tors, with the dimension four gluon condensate term illustrated in Fig. 2b. Finally, 
there are important contributions that come from the four-quark terms illustrated in 
Fig. 2c and that are the focus of this work. For light-quark mesons, the quark conden- 
sate terms are suppressed because they are multiplied by a small number, the mass of 
an u or d quark. Consequently, the most important nonperturbative contributions to 
the light-meson masses are four-quark condensates Fig. 2c, and the gluon condensate. 
Fig. 2b. 

We will denote the four-quark condensates corresponding to the correlator II'" as Q . 
These have the following form: Q = Q^'^+fl Jo|^ Q =Q''+IW,Q = (3*" + ||Jo|^ 
Q = + ^\Jo\'^ and Q = Q''*- + ||JoP, where is defined in Eq. (|l|) and 

\Jo\' ^ (0| : q7>.t-qY,qfrt% : |0). (11) 
/ 

As the four-quark condensates are mediated by gluon exchange, they contribute to the 
correlator proportional to as = 5s/4'7r ~ 0.7 at a scale of = 200 MeV Q, where 
Qs is the strength of the quark-gluon coupling. Note that we have given four-quark 
condensates corresponding to correlators for all five Dirac structures, including 11*'^, 
and that at^ refers to the current with F = ^^cr^y. Other four-quark condensates 
are needed for the correlators of the axial, tensor and axial tensor mesons, and these 
are defined below. 

As we have noted, the square of the quark condensates instead of Q appears in 
versions of the theory when the factorization approximation is applied to the four-quark 
condensates. This approximation has been successful in describing the properties of 
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hadrons at T = in free space (see, however, Ref. ||T^). Since factorization may not 
be valid when T 7^ 0, in developing our ideas we will try to avoid it as far as possible. 
Consequently, we will generally be using the exact representation of the four-quark 
condensates as given in Fig. 2c. At T = 0, the gluon condensate has the familiar value 

5,2(0| : G^G'^'^'' : |0) ^ (^^G^) ^ 0.47 GeV^ (12) 

2.2 Extension to Finite Temperature 



We want to investigate the excitation modes of hadronic matter brought into contact 
with a heat bath. Since we assume temperature equilibrium, the heat bath is the local 
environment for describing the properties of mesons propagating in the early universe 
or in the quark-gluon plasma created in the aftermath of colliding relativistic heavy 
ions in RHIC experiments. The heat bath elevates the temperature of all matter in 
the vacinity of the meson in question to temperature T, including the constituents of 
the vacuum. 

Statistical mechanics allows us to calculate the thermal average Il^{p,T), 

nr(p,T) = i [ d'xe'^-pe{xo)mj'^{x)j'^mm, (is) 



where p = (w, p) and where ((...)) stands for the Gibbs average. As in the case of the 
corresponding operator given in Eq. (^), the poles of this operator correspond to the 
excitation spectrum of the medium. Note that at finite temperature, one works with 
the retarded rather than the causal propagator. 

The quantity in Eq. ([l3|) may be evaluated in field theory using the Matsubara 



formalism |11|, which recognizes the similarity between the time evolution operator 
introduced in field theory at T = and the Grand canonical ensemble introduced at 
T ^ 0. The relevant formal distinction is the replacement of time — > i/T to obtain 
the T ^ results from the T = results. Since the limits of integration are finite and 
bounded by T at T 7^ 0, a boundary condition must be specified that is not needed in 
the case of the case of field theory at zero temperature. These differences translate into 
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the following rule for transforming a Feynman diagram at T = to a corresponding 
diagram at finite temperature: the momentum {iujn,p), where 

u!n — > 2-irTn bosons (gluons), and 

uJn 27rr(n + -) fermions (quarks), (14) 

where n runs over positive and negative integers (including zero). Thus, instead of 
integrals over w in Feynman diagrams, one has sums over n. When Feynman diagrams 
of a causal Green's function are translated in such a fashion, it can be shown that one 
recovers the Fourier transform of the advanced (retarded) real-time Green's function 
(at finite temperature T) when the analytical continuation Un — > iu>+{—)r] is made []l^ ]. 

The Matsubara formalism may be applied to the correlator in Eq. (^ and hence 
to calculate the meson masses as a function of temperature using QCD sum rules. 
One merely has to observe that the o.p.e. is made using Euclidian variables, so the 
appropriate replacement of the Matsubara frequency is now uJn ^ uj. We will be 
interested in the excitation spectrum at p=0, so we will then be making the Borel 
transform (see Sect. 2.3) on the variable uj. 

In obtaining the QCD sum rule, one must recognize that the on the right-hand side, 
the residue at the pole of the advanced (or retarded) propagator is the familiar one 
except for the presence of the occupation probability factor (for bosons) coth(mr/2T). 
On the left-hand side, the perturbative quark term Fig. 2a is evaluated in accord with 
the above considerations in terms of the Matsubara propagator 

-■ , (15) 

^^^nJo - 7 • p - m 

and the nonperturbative terms are evaluated with the temperature-dependent conden- 
sates. 

The temperature dependence of the quark condensate ((qq)) has been evaluated in 
lattice gauge theory ||T^ and in the model of Nambu and Jona-Lasinio 1 14 1 extended to 



finite temperature using the Matsubara formalism. In the latter picture, the two-quark 
condensate drops toward zero at a temperature of about 160 MeV at zero chemical 
potential. The models of the temperature dependence of the p meson mass developed in 
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Refs. [|6|, were developed using factorized four-quark condensates with ((qq)) having 
a temperature dependence of the general character found there. The temperature 
dependence or the quark condensate has also be studied in instanton models [15|. 



As we have remarked, if one avoids the factorization approximation then the actual 
contribution of ((qq)) to the light-quark mesons is small and may be ignored. One 
must instead use the four-quark condensates directly; this is the proceedure we will 
follow. As the temperature dependence of these is not generally known, we will of 
course be using the temperature dependence as it comes from our theory. 

The remaining non-perturbative quantity that needs to be specified is the gluon 
condensate. At finite temperature, we may write for the gluon condensate 

(G^) ^ (: :) + (G2)^-^- (16) 



where the perturbative term (G^)^-^- has been calculated by Kapusta |16] and found 
to be 

{G')r- - (17) 

which at Tc — 200 MeV is a very small fraction of the zero-temperature nonperturbative 
term, {G^)i^^' < G^ :). This leads to our assumption, also made in Ref. ||6|, 
that the temperature dependence of (: G^ ■.)t is weak. In the present work we assume 
that 

{G')t^{G'). (18) 

This also means that the gluon condensate is nonvanishing for temperatures T > Tc, 
so that the vacuum is still nonperturbative. 

2.3 QCD Sum Rules for Meson Masses 



In this subsection, we present the results for the QCD sum rules for meson masses Q 



extended to finite temperature using the Matsubara formalism [11|. In the next sub- 
section we will introduce similar sum rules, with however some important differences, 
for the four-quark condensates. 
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The light-quark mesons that have been succesfuhy treated in QCD sum rules include 
the lightest that couple to a given as specified in Eq. Generally, these mesons 
come in nearly degenerate isospin doublets. For the light-quark mesons corresponding 
to the elementary Dirac operators, we have the following mesons to consider: the 
fo(1370)-ao(1450) in the 0++ scalar channel (r=l); the p(776)-u;(783) @ in the l"" 



vector channel (T = 7^); ai(1260)-fi(1285) [18| in the l^"*" axial vector channel (T = 
Iblfj,)- The 7r(140)-r/(549) in the ^ psudoscalar channel (F = 75), is an exception for 
which QCD sum rules are no good. In the two scalar channels, the fg and ao have a 
different character in that only the fo couples directly to the 0"^^ scalar glueball, but 
this does not affect our calculation of the qq part of the correlator. 

The correlators 11^ (p) generally have a rich tensoral structure, and it is sufficient 
for our purposes to consider only one of these, which we will refer to as U'^ . In the 
case of the vector, axial vector, and scalar mesons n"^ is related to of Eq. (|) 
respectively by II'^II = irp'"^ g^^^ g'^^W^^ for the longitudinal correlator of the p(776)- 
a;(783); H'"^ = -lg^'''YV|% for the ai(1260)-fi(1285); and = -kW for the fo(1370)- 
ao(1450). Noting the near-degeneracy of the 1=0,1 mesons of each type except the ps, 
one can simplify the theory considerably by assuming an exact isospin degeneracy; in 
this case the correlators are multiplied by a factor of two relative to the value in the 
absence of isospin symmetry . 

We next characterize the correlators as they arise from the quark current, coming 
back to the glueball contribution of the 0"*"^ mesons later. A convenient summary of 
the various ingredients may be found in the review Q, except for the extension to finite 
temperature, which follows the arguments found above. 

The L.H.S. of the correlators arising from the quark current have the following form 
in momentum for the s, u, and av cases, 

^IhAQ\T) = arU-p') + ^r^^^ + ^r^^^. (19) 

In these expressions, —p"^ = is the Euclidian momentum, and the function In^—p^) 
characterizes the perturbative quark terms (see Fig. 2(a)) and has the following form 

/,(V) = in{-p^) (20) 
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for the s and av cases, and 



2 f°° LO 

-p^ Jo 2r 



3(-p2) J 4ml 2^ 

for the v» case. Here nsiz) = (e^ — 1)^^, np{z) = (e" + 1)^"'^, and 



). (21) 



1 /■^°^^2 
Jo -P^ + "^^ 



«n(V) = -r/ rfu;^ , , , tanh(u;/4r). (22) 



The last integral in Eq. (21) describes the damping due to thermal excitation of pions ||5[ 
and is particular to the longitudinal channel. We have assumed that the contribution 
of the continuum to the dispersion relation begins at s = Sq, and that this piece of 
the R.H.S. is transferred to the L.H.S. of the sum rule, where it regulate the large-s 



contributions of the perturbative terms. For this reason, the integrals in Eq. (22) run 
only up to 5o. This corresponds to the usual treatment of the continuum at T = 0. 

The following light-quark mesons, which have slightly more complicated currents, 
have also been successfully treated by QCD sum rules: The Di(1235) |18| in the axial 
1+- channel (J^"^'-" = ^975 9/. q)] the a2(1320)-f2(1270) H] in the tensor 2++ channel 
{JIlu = du +7i/ and the H2(1670) [|19| in the axial tensor 2""*" channel 

The correlator 11'^ for the axial meson has the same form as given in Eq. (19). For 
the tensor and axial tensor mesons, II"" has the form 

H'IhAQ^T) = arIn{-p') + br{-logi-p')){^G')+crj-^^. (23) 

The following three Q correspond to the axial (1"' ), tensor (2"'""''), and axial tensor 
(2 *") mesons respectively: Q = Q^^ — g|JoP) Q = Q^i and Q = Q"'^, where we have 
introduced rather obvious abbreviations for these channels. The appropriate In{—p'^) 
is that of Eq. (|20|) in these cases. 

The coefficients of Eqs. (19,23) are given in Table 1. We refer the reader to the 
original literature [@, 18, ^ for the details not given here. 



To obtain the sum rule, we make the Borel transform 

u'^{mIt)^Lm,u'^{Q',t)= hm ,^^^(;^)'n'r(g^r)) (24) 
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of both the R.H.S. and left-hand side. After making the Borel transform of the meson- 
mass sum rules, we may write the L.H.S. of the sum rule in the following generic 
fashion 

r 



^LhAMI.T) = ArM'^-^'\^{Ms) + Br^^^^ + ^r-^, (25) 



The function en{MB) is given in all cases by M^*-"''''^'*e„(Ms) = LMglniQ^): with 

Lu^iniQ^) ^ / duj^uj^^e ^ tanh(-^). (26) 
n! Jo 4T 

Note, in the case of the longitudinal channel for the vector meson, we have from 
Eqs. (|^) and (^ 

MleoiMB)^ du;2e «Btanh(— ) + 2 / du'nF{—)-- duj^nB{—).{27) 

Jo 4J Jq 21 S ^4^2 21 

The coefficients A, B, C, and the various n's, are given in Table 2 for all the cases. 

In the case of the 0"*"^ meson, we follow the philosophy of Ref. |jl^ in which the 
scalar meson is treated as an hybird. i.e., the current for the 0^"*" case is 

jO++=^j^ + (l_|;3|)jGB(o++)^ (28) 

where J* is the scalar quark current defined in Eq. and jGB{0++) jg ^j^g scalar 
glueball current. Details are given in Ref. |17|. Our approach differs from that of 



Ref. [P, because the mesons considered there correspond to the fo(980) and ao(980) 
mesons, which are not qq mesons. Instead, we consider the 0^^ to be the fo(1370) and 
ao(1450) mesons. Because these mesons are hybrid, there is in addition to the quark 
current and the correlator arising from the coupling of the qq pair to the glueball, 
giving rise to a "glueball" contribution 



^fGlueball ^ 4^mU — ? + 87r3(^)2r(6) + (^f^fl^. (29) 

vr TT vr 



To obtain the full correlator for the scalar meson, we combine this with the quark 
contribution in Table 2, 

n/0++ _ n2-rT'quark,0++ . i n\\2-rTlglueball 

L.H.S — P ^^L.H.S + U ~ \y\) '■^L.H.S. ' 
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where the value /3 = 0.7 reproduces the mass of the M0^"^=1370 MeV. Note that the 
cross term, proportional to — does not appear in the expression in Eq. (|^) 
because this term vanishes under the Borel transform in the form we use. 



Our correlator does not agree exactly with the one found in Ref. |17]. As in the 
case of the other mesons, we multiply the quark term by a factor of two to take into 
account that we are mixing taking the isoscalar and isovector mesons to be degenerate. 
Additionally, we use an unsubtracted dispersion relation and regulate the strong de- 
pendence on Mb in the glueball similar to the regulation of the continuum for quarks. 
Since the glueball makes a relatively small contribution for this meson, these differences 
are relatively unimportant numerically. 

To complete the sum rule for the meson masses, we need to specify the expression for 
the phenomenological correlator, ^'ii fj g ■ In all cases, H^^^ {M'^,T) is proportional 
to e~™r/^^s, where is the mass of the meson created by current . The meson- 
mass sum rule is obtained by equating the the logarithmic derivitive of the R.H.S. and 
L.H.S., 

(MlT) = (M|,r). (31) 

According to the theory, the sum rule is valid only in a region there the combination 



in Eq. (31) is relatively insensitive to Mb, which in practice means within a region 
where the sum rule has a plateau. Note that the logarithmic derivitive of the R.H.S. is 
just the square of the meson mass, and for this reason it is not necessary to know the 
constant of proportionality in I^Bjjg {M'^,T). Later, when we discuss the four-quark 
sum rules the constant of proportionality on the R.H.S. is quite important, and we will 
carefully discuss it for these cases. 

We give some of the results for mr in Table 3. We also give the values of the 
threshold parameters and the factorized values of the four-quark condensates ap- 
propriate for r = 0. As T is increased, we cannot solve the QCD sum rules for meson 
masses without knowing the temperature dependence of the four-quark condensates. 
In Sect. 3 we present two models for this. 
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2.4 Four-Quark Condensate Sum Rules 



We next want to look for a sum rule whose solution gives the four-quark condensates 
themselves. The correlators for the four-quark sum rules will be the correlators n'"(p,T) 
in Eq. (|T3| ) for the currents corresponding to the five elementary Dirac structures. 
The correlators are in most instances not the same as they were for the meson-mass 
sum rules. For example, the sum rules for 11*'^ and 11^'^ are going to be considered in this 
context. We will give our arguments in the next section why the use of the psudoscalar 
correlator the four-quark condensate may be justified while it is not justified for 171^^, 
as noted earlier. 

One can easily verify that J*"^ couples to the 1 channel from the fact that a 
nonvanishing antisymmetric tensor for the overlap in Eq. (|^ can be constructed from 
the vector meson \m{p, A)) intermediate state in ^'jiffs as 

{0\a^,\m{p,X)) = gtApA^\p)-Puel^Hp)), (32) 

where e^"^^ is the polarization vector for the vector meson. Because the current J*"^ 
couples to the same 1 channel as the p-uj, the channel is generally not considered 
as an independent case. 

To identify the quantity that we will be solving the sum rules to obtain, consider 
first the R.H.S of Il^{p) for the case of T = 0. In this case, we write for the residue at 
the pole in Eq. (^ 

^{0\J^\m^{p,X)){m^{p,X)\J^\0) = CF^{a)Pl{p) (33) 
A 

where Paip) = Ylx ^{P^ ^)'^Hp^ ^) is a projection operator whose form is specified by 
the quantum numbers for the particular meson in question, C is a overall scale com- 
mon to all mesons F, and where ^'"(q;) is a function that characterizes the differences 
among the different mesons and whose form is not needed in this work. Explicit forms 
suggested for these functions may be found in Ref. |, We give P^{p) along 

with the definitions of the correlators 11'^ (Q^ , T) for each of the five elementary Dirac 
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operators in Table 4. Note that for the vector meson, we use the transverse correlator, 

We will be making some model assumptions about the spectrum of states \moi) that 
enter Eq. , one being to replace the continuum by a set of resonances having four- 



momentum {\Ea\ = ^fP' + m'^,p). In this case, consider the matrix element {0\J^ J^\0) 
and insert a complete set of states \ma) between the two currents to obtain 

(OIJ^J^IO) = Y^J^^{0\J^\m^{p,X)){mUp,X)\J^\0) 

= Ck^ (34) 



where 



(27r)3 2Eo, ' 

where denotes contraction over the Dirac indices of P^. The value of kF is de- 
pendent on various cutoff factors that are not known, so we will eventually fix the 
value of kF phenomenologically to assure that the four-quark condensates attain their 
factorized value for T = 0. We may solve Eq. (|3^) for the scale C to obtain 



C = (0|J^J^|0)/k^. (36) 



Now, using Eqs. (jG]) and (pq ) we may write the sum in Eq. (p^ ) as 

,^ ^ (OlJ^JnO) [ d^Pa ^ pr. ,„r,^ (2vr)=^'^(p-Pa) ,,^ ^ . 
p{E,p) = J J^l^^ y'^'Po^yP) ^ d{E-Ea). (37) 



Note that with the help of Eq. ( |37D we are able to express the R.H.S. of the correlator 
n^(p) in terms of the quantity {0\J^ J^\0) that is closely related to the four-quark 
condensate ((0|qrt'^qqrt'^q|0)). Such an expression will enable us to find a set of sum 
rules for the four-quark condensates, independent of the QCD sum rules for meson 
masses. 

The extension of Eq. ( |37[ ) to the case of finite temperature is straightforward using 
the results of Sect. (2.2). We noted there that the R.H.S. for the case of finite temper- 
ature was similar to the R.H.S. at T = except for the appearance of a temperature- 
dependent occupation probability. Since the quantity appearing on the R.H.S. is a 
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meson, the correct occupation probability is coth(-^). Then, 

Pt{E,p) = -jjr^^ J -^i^l^F (a)^a(Pjcoth(— ) — (5(^-^„),(38) 

where 

In the QCD sum rule for meson masses, the higher-lying states that contribute to the 
dispersion relation (i.e., the continuum, assumed to begin at s = 5*0) are transferred 
to the L.H.S. of the sum rule, where they regulate the perturbative term at large 
momentum. We make the assumption that the continuum plays exactly the same role 
for the four-quark condensate sum rules, and therefore we characterize the continuum 
by exactly the same functions /„(M^) that characterize the meson-mass sum rules. 
Here, we parametrize Sq to be a function of temperature, the functional form carrying 
the burden of describing the dropping of the threshold Sq as the temperature is raised, 
as well as the temperature dependence of k(T) in Eq. (|39|) . With these asumptions, 
the R.H.S. of n^(p) has the following form, 

n?(..p-r)^ f W , coth(.nr/2T). (40) 

The left-hand side of 11^ {p, T) is expanded in an o.p.e. as before in the case of meson 
masses. The gluon condensate is treated as in the case of the meson-mass QCD sum 
rule. The correlator ^'lug {Q^,T) has the same form as given in Eq. (|l9|), with the 
treatment of the continuum as given in the preceeding paragraph. There is a glueball 
effect in the s channel now only in pt{E,p), since the current is purely of the form 
given in Eq. (||). The values of the coefficients are given in Table 5. 

The temperature dependence of the four-quark condensate ((Olqrf^qqrf^qlO)) is 
treated as follows. We would like to relate it to the quantity {{0\J^ J^\0)) that appears 
on the R.H.S. of the sum rule in order to obtain a set of equations to determine it. This 
connection is found by making a Fierz rearrangement of the quark operators appearing 
in ({0|qrt'^qqrt'^q|0)), so we write 

((OlqrtWf^qlO)) = 5^cr'(r)((0|jr'jr'|0)) (41) 
r' 
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making use of the following relationship among the Gell-Mann matrices 

2 

*a7*A/3 = 2(5a/3(^-^A " -^^a-y^X/S- (42) 

After making the Fierz rearrangement, the four-quark condensates Q are related to 
as shown in Table 6. We have taken | JoP = Q^: which although not an exact result 

differs from it only by terms that vanish in the factorization approximation. This is 

the only place where we make a factorization assumption (other than at T = 0). 

To obtain the four-quark sum rule, we make the Borel transform as given in Eq. (p^) 

of both the R.H.S. and L.H.S. and then equate them, 

^l.h.sXM'b^T) = n'|.^.5.(M|,r). (43) 

As in the case of the meson-mass sum rule, the equality is possible only in a region 
where the sum rule is rather insensitive to Mb- The R.H.S. and L.H.S. of the four- 
quark sum rules are given in Tables 7 and 8, respectively. The L.H.S. correlators have 
the same form as given in Eq. ( [25| ) . Note that these constitute five equations in the five 
unknowns ((0|qrt''qqrt'^q|0)) and other quantities. We will come back to the solution 
of this equation later in the next section. 

3 Toward a Theory of Four- Quark Condensates 
at Finite Temperature 

In this section, we will present two models of the four-quark condensates and the 
corresponding meson masses. For our first model, we combine the results of Sect. (2.3) 
with a parametrization of the four-quark condensates and the continuum. 

3.1 Model I: Parametrizing the Condensates and Contin- 
uum 
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We first look at a simplified model in which we use the meson mass QCD sum rules 
in the uy, av, s, a, at, and t channels, assuming that the temperature dependence of 
the four-quark condensates behaves according to 

Q^{T) = Q^{0)F{T). (44) 

For this simplified model, we will be assuming that the same function F{T) describes 
the temperature dependence of all four-quark condensates. This amounts to a factor- 
ization assumption. For the present calculation, we take 

F{T) = 1 - T^/T^, (45) 

which corresponds to the mean-field assumption of Ref. 0], where the temperature 
dependence of the quark condensate was taken to be proportional to the squar-root of 
this quantity. We take Tc = 220 MeV j^. In our second model, discussed in Sect. 3.2, 
we will be able to check factorization for T ^ Tc- 

Our second assumption regards the behavior of the temperature dependence of the 
continuum, which we are taking to begin at s = Sq{T). As we are unable to calculate 
the temperature dependence of the continuum from first principles, we parametrize 
this quantity. We have found that we get reasonable and consistent solutions of our 
sum rules using all our models if we take 

Sl{T) = Sl{0)G{T). (46) 

with 

G{T) = 1 - TVA^, with A = 331 MeV, (47) 

with S'o(O) given in Table 3. As a result of this parametrization, our results depend on 
the value of the cutoff parameter A. With this choice of Sq{T), the continuum moves 
down toward the resonance F as the temperature increases. Such behavior is expected 
to occur in a realistic system, to the extent that all the masses drop and hadronic 
resonances get increasing widths with increasing temperature. 

Using these two assumptions, and the QCD sum rule specified by Eq. (|3l|), we get 
equations for mp(T) of the following form, 

mUT) = H^{T,Sl{T),M^,Q{T)). (48) 
19 
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Figure 3: Solutions of Eq. (^81) assuming factorized forms for four-quark condensates. 

The L.H.S. of the correlators we use are given in Table 2. Note that in this model, Q 
is proportional to F(T) with the values at T = given in Table 3. The solution of 
these equations is taken in the region where mr is insensitive to M^. We have found 
that to a reasonable approximation for the av, s, and a cases this region occurs where 
attains its T = value. However, for the uy, at, and t cases, the region of stability 
in drops with increasing T. A simple, approximate description obtains by taking 
^ m,|| (T) + .5 GeV2, M^* « m„t(r), and Mj^ ^ mt{T). 
The results are shown in Fig. 3. As in Ref. Q , the mass of the vector meson drops to 
zero at T ~ Tc, but the rate of decrease is not nearly as fast here. The behavior of the 
tensor meson is particularly interesting. Between T ~ 3Tc/5 and T ~ Tc, m^{T, M^) 
changes its character several times: from a function whose relevant region of stability is 
an asymptote at large M^, to a function whose region of stability is a broad shoulder, 
and finally to a function having as its only stable solution one at relatively small values 
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of M^. We are able to express the complicated change in the location of stable regions 
by Mg « mt{T). The axial tensor meson is also an interesting case. The relevant 
region of M^* is a broad shoulder at all T. The region of stability tends to move to 
smaller , corresponding to smaller ruas, as T increases toward T = Tc- Additionally, 
in this limit the correlator becomes increasingly sensitive to the condensates. 

It is easily verified from Fig. 3 and Table 3 that the gap between the continuum 
and the meson pole decreases as expected when T approaches for all mesons except 
the tensor meson. In this case, the drop of the meson mass is more rapid than the 
dropping of Sq{T). Note that the Borel transform becomes less effective at suppressing 
features of the excited states that are not modelled in the theory as the separation 
between the resonance and the continuum decreases. Thus, in this model as well as 
the next, the theoretical uncertainties grow as Tc is approached in a fashion that is 



increasingly dependent on the unknown parameter A in Eq. (47). These uncertainties 
may be able to be decreased somewhat by introducing specific forms of {a) and 
Pql ip) given in Eq. ( |33D or by extrapolating to results obtained at asymptotically large 
T, where uncertainties may be again better controlled [p^ ]. 

Although there is considerable theoretical uncertainty in our results in the vacinity 
of T = Tc, we note the following interesting feature that may hint at new physics. 
Starting at T ~ Tc/2 the tensor meson mass begins to drop quickly, and the rapid 
drop continues as T approaches Tc. Except for this case, and that of the vector meson, 
the meson masses remain relatively large at the transition temperature. The rapid 
decrease of the tensor meson mass appears to be associated with the unusual behavior 
of the corresponding correlator at large T, noted above. Because of this, the tensor 
meson may be particularly sensitive to the temperature dependence of the four-quark 
condensate and be an important case deserving careful experimental study. 

3.2 Model II: Solving the Four-Quark Condensate Sum 
Rules 
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In this section, we remove one of the assumptions of the previous subsection, namely, 
we no longer parametrize the four-quark condensates but rather determine these con- 
densates self-consistently using the four-quark QCD sum rules whose R.H.S. and L.H.S. 
are given respectively in Tables 7 and 8. Note that these four-quark sum rules depend 
on the masses of the mesons in the av, s, ps, and t„ channels. In deriving these 
sum rules, we have already taken the meson mass on the R.H.S. to be the one corre- 
sponding to the lightest meson for each current. Thus, we take the mass my{T) to be 
that corresponding to the p — co meson; lUav (T) to be that corresponding to the oi — /i 
meson; and ms(T) to be that corresponding to the ao — /o meson. Since J*"^ couples 
to the 1 channel, we take the mass in the t„ channel to be m^(T), corresponding 
to the p — uj meson. For the mass of the pion, we take the free-space value = 140 
MeV following the observation in the model of Ref. that the pion mass is rather 
stable with temperature up to T^. The values of ms{T), my{T), and mav{T) are of 
course taken to be solutions of the corresponding three meson-mass sum rules discussed 
already in Model I, above. 

We may write the theory as follows. Note that the four-quark condensate sum rules 
may be solved algebraically for the five matrix elements ((0[ J^|0)) if Eq. ( pl) is used 
to express Q in terms of ((0| j'^jO)). The result may be written 

m,f ,f\mT) = K^{TMT)^S^{T),M'^), (49) 

i.e., as an explicit function of a Borel mass M'^ for each matrix element ((0| j'^jO)), 
the three relevant masses {m{T)), and the five threshold functions {So{T)). The re- 
maining three equations are of course the three meson- mass sum rules of Eq. (^) , 

- r 

ml{T) = H^{T,mr,S^{T),M^,Q (T)) for E = s, v, and av. (50) 

To get a completely self-consistent set of equations, we again express Q in terms of 
((OIJ^J^IO)) using Eq. (|l|). 

Our Model II thus consists of eight coupled nonlinear equations for the eight un- 
knowns, consisting of the five matrix elements ((0[ J^|0)) and the three meson masses 
nis (T) , rriv (T) , and niav {T) . As we are unable to calculate the temperature dependence 
of the continuum from first principles, we continue to parametrize this as in Eq. (^) 
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Figure 4: Temperature-dependence of R(T) defined in Eq. (^) (solid line) for the scalar 
four-quark condensate compared to temperature dependence of four-quark condensates in 
the factorized, mean- field approximation (dashed line). 



and Eq. (^), above. We take the thresholds Sq to be the same as those in the T = 
meson-mass sum rules, Table 3, except in the case of the ps and to- channels, where we 
take Sq" = and S^' = 1300 MeV. The parameters in Eq. ( p9| ) are obtained by 
solving the four-quark condensate sum rules at T = and adjusting the values of these 
parameters to obtain the known factorized values of the matrix elements (0| J^J^|0). 
Note that we are assuring that the quark condensates attain their free-space values by 
our choice of kF . The values that assure this are given in Table 9. 



To obtain numerical solutions of the nonlinear system given in Eqs. (49, 5C), we 
begin at T = 0, where we know that the solution of our equations gives the well-known 
results for the meson masses and (factorized) four-quark condensates. We then slowly 
increase the temperature, letting T — > 5, solving by iteration the eight equations in 
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T/L 

Figure 5: T-dependence of QCD sum rule determined four-quark condensates as a function 
of temperature. Note that the T-dependence is given relative to the the scalar four-quark 
condensate shown in Fig. 4. 

eight unknowns until a stable solution is reached for the for the new meson masses and 
condensates at T = 5. Then, we increase the temperature T + 5 and iterate again 
for a stable solution. The process is continued up to the transition temperature Tc. 
Our results do not rule out the existence of multiple solutions at finite temperature. 
We have in all cases presented the solutions that evolve continuously from the solutions 
at T = 0. 

We obtain the solution of each sum rule by finding the region of stability as a 
function of the Borel mass for each quantity of interest. We thus find a Borel mass for 
each of the eight unknown quatities at each T. Once we know the condensates, we can 
solve the sum rules in the a, t, and at channels. Clearly, the meson masses in these 
three channels are determined as a consequence of the eight coupled equations, but the 
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Figure 6: Meson masses as a function of temperature using our QCD sum rule determination 
of the four-quark condensates. 



corresponding mesons do not appear as part of the four-quark condensates sum rule. 

We give the numerical solutions for the coupled equations next. We will denote the 
the temperature dependence of the condensates by R^{T), so that 

((oiJ^j^lo)) = i?^(r)(o|j^j^|o). (51) 

First, in Fig. 4, we show the behavior R^{T) for the four-quark scalar condensate. 
The solid curve shows the solution to our coupled equations, and the dashed line shows 
for comparison the mean-field temperature dependence given in Eq. (^). We see that 
the scalar four-quark condensate falls off with temperature very much as expected 
from factorized approximations 0. We have found that the details of the temperature 
dependence in Fig. 4 is governed by our choice of G{T) in Eq. (|47|). 

Next, consider the relative temperature dependence of all of the four-quark con- 
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densates. These are shown in Fig. 5. We see here that the scalar condensate has the 
slowest fall-off with T. Additionally, it is clear that all condensates T fall off in T 
with different rates; the ratio R^{T)/R'^{T) is much less sensitive to G{T) than R^{T) 
itself. Vacuum saturation would imply the same rate of fall-off, so we see that even at 
a modest temperature of T ~ Tc/2 there is a considerable violation of factorization. 

At r=0 the factorization approximation of the four-quark condensates is also in 
question. In a study of the ratio of the isovector hadronic to the muon pair production 



in e'^e~ data [10| using a sum rule that weights the heavier states more than the usual 
sum rule a value of the vector four-quark condensate is larger than the factorized 
value. Using this value we find that we get a satisfactory fit to the rho-meson mass. 
Moreover, in the present work, we observe that the T = 0, B = sum rules give 
all meson masses consistent with experiment for our parameterization with T = 
factorization. We have explored the additional features of four-quark condensates for 
finite T, and there is no difficulty for us to modify our results to include correction 
factors for T = 0. 

We give the numerical solutions for the meson masses in Fig. 6. Note that the 
temperature dependence of the meson masses for all cases except the tensor and axial 
tensor mesons are quite similar to those of Model I shown in Fig. 3. Because the tensor 
and axial tensor show more sensitivity to the four-quark condensates, these cases may 
provide particularly interesting candidates to explore experimentally. 

We have noted that accuracy of the correlator in the psudoscalar channel for the 
mass of the tt meson is very poor. Novikov, et al. |24| have indicated that for channels 
that couple easily to the vacuum, such as the psudoscalar channel, instanton correc- 
tions identified as "direct instantons" may be particularly important. The subject of 
the instanton contributions to meson mass correlators has studied in some detail in a 
recent review [p^]. Here different models of correlators have been compared, showing 
how instantons can improve the correlator out to large separations r. Although direct 
instantons presumably also contribute to the psudoscalar and other four-quark conden- 
sate sum rules that we have used, we have normalized these sum rules to the factorized 
value of the four-quark condensate at T = 0. The fact that we are only relying on 
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the four-quark sum rules to provide the temperature dependence of the condensates 
and not their T = values presumably improves the reliability of the sum rule for our 
purpose. 

4 Summary and Conclusions 

This paper has focused on the physics of the four-quark condensates Q^, which 
constitute fundamental, nonperturbative features of quantum chromo dynamics needed 
for understanding how hadron masses (especially those of the light-quark mesons), as 
well as interactions among hadrons, arise. These as well as the other condensates all 
govern essential details of the phase transition that is believed to have occured as the 
early universe expanded and cooled from the era dominated by the quark-gluon plasma 
into the next, in which the quarks and gluons had coalesced into the familiar hadrons. 
Understanding the physics of this phase transition has been identified as one of the 
opportunities that presents itself to nuclear and particle physics community at the next 
generation of experimental facilities of RHIC and perhaps the LHC. 

In spite of the importance oiQ^ = ((Olqrf^qqrf^qlO)), these quantities are currently 
quite poorly understood at both the phenomenological and at the theoretical levels. 
For example in the literature, the temperature dependence of is almost universally 
described through a factorized, mean-field approximation in which cx ((qq))'^ oc 
1 — T^/T^, with Tc ~ 200 MeV. In such a formulation, the dominant nonperturbative 
quark contribution to the masses of light-quark mesons is the quark condensate ((qq)). 
In this paper, we have proposed a dynamical theory for the temperature dependence 
of that generalizes the familiar gap equation for the quark condensate ((qq)). In 
our theory, the masses of the light-quark mesons are practically independent of ((qq)); 
instead these masses depend on four-quark condensates that do not factorize for T ^ 0. 

Our work is based on finite-temperature QCD sum rules in the Matsubara formalism 
and requires a model of the hadronic spectrum. The model of the hadronic spectrum 
we have chosen is the discretized model adopted in Ref. |^. It contains, in addition, a 
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parameter that characterizes the temperature-dependence of the continuum. Equations 
are obtained couphng light-quark meson masses mr and ((0|qrt^qqrt^q|0)). Stable 
solutions are found in which ((0|qrt^qqrt^q|0)) all decrease as T increases toward Tc, 
but at different rates, so that factorization becomes increasingly violated the farther T 
is varied away from T = 0. 

We stress the connection between the temperature dependence of meson masses 
and the fundamental four-quark condensates in hot quark matter. The tensor 2++ 
mesons a2(1320)-f2(1270) are particularly important for giving quantative information 
on the four-quark condensates. The temperature dependence of meson masses consti- 
tute predictions of our theory that may at least in principle be checked experimentally 
by measurement of final electromagnetic states (77, e"'"e~, or or perhaps semi- 

electromagnetic states such as n^. 

Our results do not depend on achieving thermodynamic equilibrium in rclativistic 
heavy ion collisions, but only local thermodynamic equilibrium. We also note that the 
temperature- (and more generally, density-) dependence of meson masses may affect 
predictions of transport models modelling the evolution of matter into and out of the 
quark-gluon plasma. Results such as the ones we give may be incorporated in such 
models able to follow the evolution of the local temperature. 
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Table 1: Meson mass correlator ^'[^n.s.iQ'^^'^)- Coefficients correspond to Eq. (19) for the 
av, s, and f y cases, and to Eq. (23) for the a, t, and at cases. For the s case, there is also a 
glueball contribution, given in the text. 



Case r n or n(r, b) br n{T, c) cr 



av 


1 


+ 


1 


TT 

12 


2 


277^0;, 


s 


1 




1 


TT 
8 


2 




V\\ 







2 


TT 

12 


3 


-27r2Q;, 


a 


1 


— (l + — ) 

StT ^ TT ) 


1 


TT 

24 


2 




t 


2 


StT ("^ TT ) 




Stt 
9 


1 




at 


2 


StT ("^ TT ) 




Stt 
9 


1 
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Table 2: Borel transform of meson mass correlator U.'[ jj ^{Mq^T). Coefficients correspond 
to Eq. (25) in all cases. 



Caser n Ar n{r,B) Br n{r,C) Cr 



av 


1 


_i_ 


+ 





TT 

12 


1 


271^0;, 


s 


1 


- ( 







TT 
8 


1 


TT^as 







J_ 




1 


TT 

12 


2 




a 


1 


1 

Stt 


(1 + ^) 





TT 

24 


1 




t 


2 


_3_ 
Stt 


(i-f) 


-1 


Stt 
9 







at 


2 


_3_ 
57r 


(i-t) 


-1 


Stt 
9 
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Table 3: Results for meson-mass sum rules, T — 0. Factorized values of the four-quark 
condensates are also given. 



Case r 5*0^ (MeV^) rm (MeV) <Q >fact / < qq >' 



av 1700 1140 



3000 1360 



1500 748 ^ 



2740 1240 



2500 1270 -f 



at 3500 1580 f 
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Table 4: Definition of projection operator P^{p) and four-quark condensate correlator 

n'r(g^T). 



Case r P^ip) Def. n'^ 



av 



s 1 7rn^ 

- PuPpQuX - PvPxQup) 

ps -1 nUP' 



Table 5: Four-quark condensate correlator n^ji^5 ((5^, T). Coefficients corresond to Eq. (19). 



Case r n ar n(r, b) br n{T, c) cr 



s 



av 2 + - 1 ^^a 

1 1 f 2 leTT^a, 

1 + 1 -f^ 2 27r'a, 
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Table 6: Fierz rearranged four- quark condensates. 



Case r 


- r 

Q 






Cs 






V 




29 
27 


5 
9 


14 
9 


14 
9 





av 


Q'^ + IW 


11 
27 


11 
9 


22 
9 


22 
9 





s 




2 
9 


2 
3 


22 
3 


14 
3 


1 

3 


U 




53 
108 


19 
36 


11 
9 


4 
9 


1 
4 


a 




25 
54 


11 
18 


5 

18 


25 
18 


1 
4 


ps 


Q«* + i|Jo|' 


2 
9 


2 
3 


22 
3 


14 
3 


1 
3 


t 




1 
3 


1 


-2 


-2 





at 




1 


1 
3 


2 


2 
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Table 7: Four-quark correlator n^^^ (M|,r). 

mJ-jy)) rr^^ ^^^^ 

s mj^j;m ^ ^^th (i^) e— ^/^i 

_ ((o|jy;-|o)) ^ ^oth (^) e-?/^I 
- <<°l-^JI°» 7r coth (^) e-^l^l 



Table 8: Four-quark correlator ^l,h.s.{^b^T). 

CaseP n Ar n(r, 5) Br n{T,C) Cr 

av 2 + - ^TT^a 



I f 1 leTT^a, 



s 



1 ^(l + f) -fL 1 27r2a 
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Table 9: Values of k,^ and factorzied values of < 0\qTqqTq\0 >. 



Case r (GeV^) (0| : qVqqVq : |0)/(0| : qq : |0)2 
av 0.0021 I 
s 0.00073 ^ 
U 0.028 -1 

ps -0.00034 ^ 
vi 0.0036 
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